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I. INTRODUCTION 

This paper describes numerical techniques that improve the accuracy and stability of 
algorithms for the small-disturbance and full-potential equations used to calculate 
transonic flows. For the small-disturbance equation, the algorithm improvements are: 
1) the use of monotone switches in the type-dependent finite-differencing, and 2) the 
use of stable and simple second-order-accurate spatial differencing. These improve- 
ments are for steady and unsteady transonic flows. For the steady, full-potential 
equation, the improvement is in the use of a monotone switch in the type-dependent 
finite-differencing of an approximate-factorization (AF2) algorithm. All these 
improvements can be implemented in present computer codes by making minor coding 
modifications. 


II. SMALL-DISTURBANCE EQUATION IMPROVEMENTS 
1 • Monotone Implicit Algorithms 

Most of the algorithms in transonic small-disturbance potential codes use the Murman- 
Cole [I] method of switching the differencing scheme for regions of supersonic and 
subsonic flow. A weakness in this method is that it allows stable solutions for flow 
fields containing entropy-violating expansion shocks in one-dimensional model cases, 
as shown in Ref. 2. In comparison, the implicit method presented here, which uti- 
lizes the monotone Godunov switch [3], does not allow such nonphysical solutions. In 
two-dimensional flows over airfoils [2], this weakness in the Murman-Cole switch allows 
numerical instabilities to develop in the calculations near the leading edges of the 
airfoils. In some cases of steady flows, converged solutions are obtained by the 
monotone approximate-factorization algorithm (MAF-G) (see Fig. 1 from Ref. 2), whereas 
calculations using the Murman-Cole switch (AF2) are unstable. In some cases of 
unsteady flows , the monotone, method (LTRAN2-MG) allows the use of larger time-steps 
than the older method (LTRAN2), which uses the Murman-Cole switch. In Figs. 2(a) 
and 2(b) [2], a case is shown in which both methods use the larger time-step and the 
older method diverges (i.e., computer overflow); here the older method required a 
time-step 12 times smaller in order to agree with the monotone method. The jaggedness 
near the leading edge in the LTRAN2-MG results is due to the use of measured ordinates 
from an experimental model. This jaggedness provides a severe test case for the 
stability of any algorithm. 

Two-Dimensional Steady Algorithm 


Governing Equation 

The MAF-G algorithm is for the low-frequency, unsteady, small-disturbance potential 
equation [2] 

2kM*<j> = [(1 - M 2 ) - (y + l)rf"$ ]<j> + <J> (1) 

°° T xt °° «> x xx yy 

Although Eq. (1) is physically meaningful for unsteady flows, MAF-G is nonconservative 
in time for more efficient calculations of steady flows. The algorithm for unsteady 
flows, namely, LTRAN2-MG, is presented in detail in Ref. 2; it was used to generate 
the results shown in Fig. 2. 


Algorithm 


MAF-G is a modification of an approximate-factorization implicit algorithm, AF2 , which 
uses the Murman-Cole switch. The MAF-G scheme uses the monotone switch of Godunov. 

It is given by the following two-step, finite-difference approximation to Eq. (1) at 
mesh point (i,j). 
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and where o> is a relaxation parameter, R n is the residual given by 


R li - 'vM + 




and u is the sonic value of <J> X . The Murman-Cole switch only uses the Ei,j switch 
[2]. The MAF-G scheme is identical to AF2 except in regions where the flow field 
changes type — near the sonic lines and shock waves. 


2. Second-Order-Accurate Supersonic Spatial Differencing 

Current methods for calculating transonic flows with the small-disturbance equation 
typically are only first-order accurate in the supersonic regions of the flow. How- 
ever, calculations using the full-potential show significant improvements in accuracy 
when second-order methods are used. In this paper, a stable, simple algorithm [4] 
is described that is a second-order-accurate extension of the implicit monotone 
algorithm described above. For steady flow. Figs. 3(a)-3(c) show calculations of 



flow over a Korn airfoil, for coarse, medium, and fine grids. Note in Fig. 3(a) the 
improvement in resolution with the second-order method. Also, the convergence rates 
of the two methods are essentially the same. The improvement in the unsteady algo- 
rithm (LTRAN2-MG) and improvements in unsteady flow calculations are presented in 
Ref . 4 . 


Algorithm 

The modifications to MAF-G to implement second-order accuracy are made by changing 
only the first step of MAF-G, given by Eq. (2). Let 
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Then the modified step 1 is given by the following: 

Step 1': 
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Now the residual is given by 
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III. FULL-POTENTIAL EQUATION IMPROVEMENT 


Monotone Implicit Algorithm 


Most of the algorithms in transonic, full-potential codes use type-dependent differ- 
encing that is a generalization of the Murman-Cole switch, such as the method of 
Jameson and the AF2 method of Holst and Ballhaus [5]. These methods suffer from a 
weakness similar to that of the Murman-Cole switch in dealing with nonphysical expan 
sions shocks. A monotone implicit method (AF2) has been developed [2] that elimi- 
nates this deficiency at sonic expansion points, as demonstrated by the calculations 
shown in Figs. 4(a) and 4(b) [2] for a two-dimensional case of flow over a modified, 
double-wedge profile. The full details of the MAG algorithm are given in Ref. 2. 
Here we describe the essence of the new switch. The AF2 method uses a nonmonotone 
switch, called' upwind-density biasing, to stabilize the calculations in supersonic 
regions of the flow. That switch is implemented in the mass flux. An example, for 
P<P X , is 


P i+(l/2),j 6 X <f, i,j ~ {P i+(l/2),j ' V i,j tP i+(l/2),j P i-(l/2),j ]} Vi,j 

MAF employs a monotone swLtch based on Godunov's ideas. An example, for p<j> x , 

p i+(i/ 2 ) ,j 6 x^i, j ^ P i+(i/ 2 ),j \+(i/ 2 ) , j '■ P i+(i/ 2 ) , j p *^^x < * > i,j 

+ { V(x/2).;jl p i-(l/2).J - P * ]} Vi,j 


(5) 

is 


( 6 ) 


The essence of the difference between MAF and AF2 is seen by comparing Eqs. (5) 
and (6). Whereas Eq. (5) uses the switch v to smoothly shift the flow variables 
upwind as the flow becomes supersonic, Eq. (6) uses the switches v and v to 
smoothly shift the flow variables to sonic values p* and q*, and to shift flow 
variables located upwind away from sonic values. Both methods use rotated differ- 
encing and treat shock waves identically [2], but the monotone method is more stable 
at sonic expansion points. 
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Fig. 1 


Converged solution obtained from algo- 
rithm using the monotone Godunov 
switch MAF-G; algorithm using nonmono- 
tones switch AF2 diverged. 
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a) kt = 75° 



Fig. 2 Comparison of algorithms using the monotone Godunov switch LTRAN2-MG and 
nonmonotone switch LTRAN2. Plots of upper-surface pressure coefficients 
of a NACA 64A010 airfoil (experimental model ordinates) in pitching 
motion. Time given in degrees of motion kt and Moo = 0.80. 
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a) Coarse grid b) Median grid c) Fine grid 

Fig. 3 Comparison of first- and second-order methods for steady flow over a Korn 
airfoil; pressure-coefficients plot: Moo = 0.755 and a Q = 0. 


MODIFIED DOUBLE WEDGE PROFILE 




a) Converged solution using MAF, b) Iteration No. 21. 
nonmonotone algorithm diverged. 

Fig. 4 Comparison of monotone MAF and nonmonotone AF2 algorithms for flow over 
a modified, double-wedge profile. 
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